Landau levels in asymmetric graphene trilayers 
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The electronic spectrum of three coupled graphene layers (graphene trilayers) is investigated in 
the presence of an external magnetic field. We obtain analytical expressions for the Landau level 
spectrum for both the ABA and ABC - type of stacking, which exhibit very different dependence 
on the magnetic field. The effect of layer asymmetry and of external gate voltages can strongly 
infiuence the properties of the system. 
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I. INTRODUCTION 



The extraordinary level of interest in the study of single 
layer graphene has led to the prediction and observation 
of several unusual phenomena not found in other low- 
dimensional systems. These new properties are mainly a 
consequence of the chiral and massless character of the 
quasiparticles in graphene^. This new material is not 
only expected to lead to several technological applica- 
tions, but has also helped shed light on relativistic quan- 
tum effects, such as Klein tunneling^''^ and zitterbewe- 
gung. There is currently a search for experimental meth- 
ods for producing high quality samples of graphene in 
large quantities. However, the current experimental tech- 
niques can also create carbon structures with two or more 
layers. It has already been recognized that graphene bi- 
layers (i.e. two coupled layers of graphene) can display 
interesting new properties that are distinct from those 
of single-layers^ and that can also be eventually har- 
nessed for the development of devices. One important 
aspect of these structures is the fact that, in comparison 
with the in-plane interactions, the comparatively weaker 
interlayer coupling can yet exert a significant influence 
on the carrier spectrum. Thus, whereas the electronic 
dispersion at the vicinity of the Fermi energy in single 
layer graphene is linear, in bilayer graphene it displays 
an approximately parabolic shape with the appearance of 
higher energy bands. Moreover, in striking contrast with 
single layers, the electronic spectrum of bilayer graphene 
has been shown to develop a gap in the presence of an 
external electric field. Thus, the interlayer coupling in 
stacked layers of graphene gives rise to a rich set of prop- 
erties that are not found in monolayers, and can be ex- 
pected to be of particular significance in structures with 
three or more layers. 

In this work we investigate the properties of three cou- 
pled layers of graphene, i.e. trilayer graphene (TLG) in 
the presence of an external magnetic field perpendicular 
to the plane of the layers. The properties of TLG in the 
absence of a magnetic field have been considered in the 
literature within a tight-binding model (see, e.g.^'^), as 
well as through first principles calculations®. The effect 
of an external magnetic field was calculated by means 
of an approximation based on the mapping of stacked 



graphene layers to an ID tight-binding chain by Guinea 
et al.^ . Recent experimental studies have investigated 
the Landau level spectrum^*^'^^ and the magnetoconduc- 
tance of TLG^^"^**. These results showed that one impor- 
tant aspect of TLG is the fact that the energy bands at 
the vicinity of the Fermi energy are very sensitive to the 
particular type of stacking of the layers. The two more 
relevant stacking rules are the rhombohedral, or ABC 
stacking, and the Bernal, or ABA stacking. In each case, 
the relative positions of the lowermost layer (C in one 
case and A in the other) helps dictate the possible sym- 
metries of the subsequent wavefunctions associated with 
each layer. The type of stacking is also relevant to the 
properties of the Landau levels of the TLG, and it has 
recently been reported^^ that 15% of exfoliated TLG has 
rhombohedral (ABC) stacking. The goal of the present 
paper is to present analytical results for the spectrum 
of TLG in a magnetic field considering different poten- 
tials in each layer. In order to do that we perform a direct 
diagonalization of the six-band continuum model and ob- 
tain analytical expressions for the Landau level spectrum 
as function of magnetic field and the potentials at each 
layer, for both the ABC and ABA stackings. In partic- 
ular, we calculate the TLG spectrum in the presence of 
electric fields that break the layer symmetry. It has re- 
cently been shown that these different potentials can lead 
to the opening of a gap in the TLG spectrum^'^^"^®. In 
the present paper we focus on the modifications of the 
Landau levels as function of the layer potential. 

The paper is organized as follows: in section II we 
present the model and solve the resulting system of equa- 
tions for the ABA stacking, followed by section III in 
which we present the model and obtain solutions for the 
ABC stacking. Section IV shows and discusses the cal- 
culated Landau level spectra and finally, the results are 
summarized in the conclusions. 



II. ABA STACKING 

Let us consider a system consisting of three cou- 
pled graphene layers, in the context of the tight-binding 
model. We assume nearest neighbor hopping between 
sites within each layer, described by the coupling pa- 
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FIG. 1: Diagramatic scheme of couplings in graphene trilayers 
for ABC (a) and ABA (b) stackings. 



rameter ig- In the continuum approximation, this pa- 
rameter determines the magnitude of the Fermi velocity 
vp = tQah~^^/3/2 fa 10^ m/s. The nearest neighbor in- 
terlayer coupling scheme is Ai — B2 — A^, with coupling 
parameter ti (see Fig. 1). In order to capture some of 
the finer details of the spectrum, we also include a remote 
coupling term ^2 between sites of sublattices Ai and A^ 
and is between Bi and B3. The Hamiltonian is given as 
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where it ^ px + ipy, with p^^y being the components 
of the in-plane momentum; C/1,2,3 is the potential in 
each layer, respectively, Uq is the onsite energy at sub- 
lattices Ai, B2 and A3, and we defined the eigenstates 
as^=[V'Ai, iipBi, ipB2, iipA2, i>A3, iV'ssJ^-In 
the presence of a uniform magnetic field in the z direction, 
with the gauge A = (0, Bx, 0) and, for a given sublattice 
L, tpL{y) = 
equations: 



^gjfeyy^ one obtains the following system of 



A^(j)Bl + t'i<l>B2 + t'2(f>A3 = (e - Ul - Uo)4>Al, (2a) 
A~(j>Al+t'3(f>B3 = ~'if'-Ui)(l)Bl, (2b) 

A'~(j>A2 + t[(j>Al + t'i(j)A3 = (e - M2 - Mo)</'S2, (2c) 
A+(j)B2^-{e-U2)(j)A2, (2d) 

A^(j)B3 + t'i(f>B2 + t'2<t>Al == (e - U3 - Uo)0A3, (2e) 
^~0A3 + i3'^Bl = -(e-M3)0S3, (2f) 



where e — E/hvp, Ui = Ui/hvp, t[ = ti/hvp and 
eB/hvp and we defined the operators 



A^ 



dx 



± {ky - l3x), 



(3) 



which obey the commutation relation [^+,^^] = 2/3. 

For Ui = U2 = U3 = U the system can be easily 
solved by making use of its reflection symmetry. Thus, we 
can define symmetric and antisymmetric combinations of 
the spinor components. For the antisymmetric case we 
obtain 0g = ■^{<t>Ai - ^As), and 0// ee -^(^si - </'B3)- 
That leads to the following pair of coupled equations 



A <I>G = -{e-u-t'^] 



IJH, 



A^(i)H = {e-u + t'2~ ua)(f)c 



(4a) 
(4b) 



For the sake of convenience, let us now define the opera- 
tor 



Z = A-A+ = ^-ih- Pxf - p. 
We can now decouple the equations to obtain 



(5) 



{Z + 2p)4>G^-[{e-u'f-{5uf]<j)G, (6) 

which corresponds to the equation that gives the spec- 
trum for a single graphene layer under an effective elec- 
trostatic potential U' = C/ -f (ts — t2 + Uq)/2 as well 
as a finite gap term given by 5U — {Uq — t2 — ^3)72. 
Thus, the Landau level spectrum in this case is given by 
e == u' ± y/2(3n + {6u)^ . From Ref. [15] we have ^2 = 0.04 
eV, ts = -0.02 eV and Uq = 0.05 cV. Hence, the addi- 
tional effective potential leads to an energy shift of —5.0 
meV, and SU — 15 meV. Thus, the effect of the remote 
coupling terms t2 and ts, as well as the on-site energy 
term Uq is only to introduce a small shift of the disper- 
sion branches and to generate a small gap in the energy 
spectrum. 

For the symmetric case we have 4>c = -y^{4>Ai + (t>A3) 
and <j)n = ^(<^bi + 0S3)- The equations become 



A^ 



A- 



OD 



V2t[(j)B2 = {e-u-t'2-uo)(jJc, (7a) 



= -(£-^+4; 



(7b) 



and 



A (l)A2 + V2t[(l)c ^ ie-u-uo)(j)B2, (7c) 



A^ 



'JB2 



= -{e-u) 



»A2, 



(7d) 



These equations can be decoupled, resulting in the 
fourth-order differential equation 



{Z^ + \lZ-\2}^C=Q. 



(8) 



where Ai = (e — u — uo)(e ^ u) + [e — u + t'^){e — u — uq - 
tj) -1-2/?, and A2 = ~~{e~u — UQ){e — u){e — u + t'^){€ — u- 



Uo — t'2)+2{t[)'^{e — u){e — u + t'.^). This equation is similar 
to the one describing bilayer graphene. A second-order 
equation can be obtained by calculating the roots of the 
second-order equation as 



^2 



{Z-z+}{Z-z_}<j>c=0, 
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z± 



Ai 



± 



Ain2 



A, 
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(10) 



where we set ky ^ 0, since this term only introduces a 
shift of the wavefunction In particular, for uq, ^2.3 = 
the equations yield results that are identical to those of a 
gapless single-layer and bilayer graphene, i.e. the Landau 
levels are found as the solutions of 2/3(n + 1) = z±. As in 
the previous case, the addition of remote coupling terms 
introduces a gap in the spectrum. 

A more realistic description of TLG structures should 
take into account asymmetries between the different lay- 
ers, which can be brought about by the interaction with 
a substrate or by gating. In order to assess the effect 
of layer symmetry breaking in the spectrum, let us now 
consider the case Ui ^ U2 ^ U3. In addition, we now 
consider t2 = t^ = Uq = 0, since we assume that the 
shifts caused by the layer potentials are more significant 
than the effect of these terms. A simple substitution al- 
lows us to write 



[A+A- + {e- uif](l)Ai = t'i(e - "1) 



»S2, 



[A+A- + (e - U3f](l)A3 = t[(e - U3)(I)B2- (11) 

In addition, we also have 

[A-A+ + (e - U2f](f>B2 = t'i(e - "2)(0A1 + 0A3). (12) 

As in the previous case, we introduce symmetric and anti- 
symmetric combinations of wavefunctions, and let us also 
define A = (ui — U3)/2, s ^ (ui + u:i)/2 and 5j ~ e — Uj, 
j — 1,2,3, in order to simplify the notation. Thus, af- 
ter some algebra, we can obtain the following 6th-order 
differential equation for (J)b2 as 

{[A-A+ + dl + 2p][A-A+ + S2^][A-A+ + S3^ + 2p] 
-tfS2{e-s)[A-A++di^ + 2p] 
-tfS2{e - s)[A-A+ + Ss^ + 2p] 
+4tfAH2{e-s)'j^B2^0. (13) 

It is seen that for Ui = U2 = U3 (i.e. Si = 62 = S3, A = 
0), we recover the previous solutions. One can rewrite 
Eq. (22) as 

[Z^ + aiZ^ + a2Z + a3]0B2 = 0. 



with the Z operator defined above and 



a\ 



Si + Si + Si + 4/3, 



(14) 



(15a) 



{S\ + 2ml + 2/?) + {S\ + 2/?)j2 + [Si + 2/3)(5i 



-tf(52((5l+(53), (15b) 

^3 = {Si + 2p)Sl{Sl + 2/3) - 2/3if ,52(<5i + -53) 

-tf<5i(52(53('5i+(53), (15c) 

This equation can be written as 

{Z-Zi}{Z-Z2}{Z-Z3}(/.B2=0, (16) 

where Zj, j = 1,2,3 are the three roots of the cubic 
equation, Eq. (14). Therefore, the spinor component 
(l)B2 is a solution of 



9S2 



dx^ 
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fixfcj)B2 



-{z,+m 



B2- 



(17) 



For zero magnetic field, this equation allows us to obtain 
plane wave solutions for each dispersion branch. The dis- 
persion relation can be obtained by setting Zj = —k"^. It 
can be immediately seen that the energy gap at fc = 
can be found by solving the equation a3 = 0. For fi- 
nite magnetic fields, the solutions are expressed in terms 
of Hermite polynomials. Therefore, for the Landau lev- 
els we obtain the relation Zj = —2j3{n + 1). Thus, the 
energies are found by solving the algebraic equation 

-Win + l)f + ai [2/3(n + 1)]^ - a2[2l3{n + 1)] + ^3 = 0. 

(18) 
It is evident that for Ui = U2 = U^, we have A = and 
the last term of Eq. (13) vanishes. The spectrum should 
then consist of a superposition of the spectra of single 
layer graphene and bilayer graphene. 



III. ABC STACKING 

Let us consider three coupled graphene layers in the 
ABC stacking configuration. For the sake of simplicity, 
let us retain only the nearest-neighbor coupling terms. 
In this case, the Hamiltonian can be written as 
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where Ui^2,3 is the potential in each layer, respectively, 
and we defined the eigenstates as before. 

Thus, one can obtain the following system of equations: 
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(20a) 
(20b) 



A-(j)A2 + t'(t)Al = (e - U2)(t)B2, (20c) 

A+4>B2 - t'(f)B3 = -(e - U2)<^A2, (20d) 

A+<l>B3 = {e-U3)(l)A3, (20e) 

A-^A3-t'(PA2 = -ie-U3)4)B3, (20f) 

where e = E/hvp, ui — Ui/hvp, t' — t/hvp and l3 — 
eB/hvp- In order to decouple these equations, let us 
first obtain (f)Bi and (pAz in terms of (pAi and (pB^ from 
the second and fifth equations as 
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(21) 

and substitute these expressions in the first and sixth 
equations, respectively, to give 

A+A-(t>Ai - t'{e - ui)(j)B2 = -(e - uif(l)Ai, (22a) 

A-A+^B3 - t'{e - uz)(j)A2 = -(e - u:if(j)B3. (22b) 

Equations (22a) and (22b) allow us to obtain (t)B2 and 
0A2 in terms of 4>ai and 0b3, respectively. Thus, by 
substituting them in Eqs. (27) and (28), respectively, 
and after some tedious algebra, one can obtain a 6th 
order differential equation as 

I [A-A+ + 5l + 2p][A-A+ + 5l][A-A+ + 5l - 2p] 
-t'^625:,[A-A+ +5l + 2p] 

^t'^5iS2[A-A+ +5l- 2p] + t'^5i5:,^(l)Ai - 0. (23) 

It is interesting to compare Eqs. (13) and (23). The 
former remains invariant if one switches the potentials 
in layers 1 and 3. Equation (23), on the other hand, is 
found to be invariant under an interchange of potentials 
between the topmost and lowest layers together with a 
reversal of the magnetic field. This reflects the different 
symmetries of each stacking of TLG. 

As before, we can obtain the Landau level spectrum 
by rewritting Eq. (23) as 



with Z = A^A+, and 



(24) 



71= 61 + 51 + 51 (25a) 

72 = {51 + 2P){51 - 2/3) + {5! + 2/3)51 + (Sl - m^l 

-{t'f{5i+ 5^)52, (25b) 

73 ^ -{t'f52M5l + 2(3) + 5,{5l - 2(3)] + {t')%5, 

+ (Si5253)^ + 2P{51 - 5f)5l - A{(352)^ (25c) 



This equation can be written as 

{Z - Zi}{Z - Z2}{Z - Z3}<()A1 



0, 



(26) 



where Zj, j = 1,2,3 arc the three roots of the cubic 
equation, Eq. (24). Therefore the spinor component (f>Ai 
is found as a solution of 



'JAl 



dx^ 



+ {ky - (3x)'^(/3Ai = -{Zj + /3)(/)^i 



(27) 



For the particular case oi Ui — U2 ~ U3 ^ and zero 
magnetic field, we can obtain plane wave solutions by 
setting Zj — — fc^, where k is the in-plane wavevector. 
Thus, Eq. (24) can be rewritten as 

e^ - (3fc2 + 2t'2)e4 + (3fc4 _ 2kH'^ + 1'^ - fc^ ^ 0. (28) 

Let us now consider the low-energy limit e << t' . That 
allows us to neglect the higher-order powers of e to obtain 



fc3 
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t'2 ^1 - 2fc2/i'2 + 3fc4/i" 



(29) 



Thus, for small wavevectors, the dispersion relation in- 
creases with the third power of k. 

The Landau levels can be obtained using the relation 
Zj ~ —2(3{n + l)^ which leads to the algebraic equation 

-[2/3(n + 1)]3 + 71 [2(3{n + l)f - 72[2/?(n + 1)] + 73 = 0. 

(30) 
As seen above, for zero potential in each layer, we can 
find a simpler algebraic relation for the energy, namely 
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+ [l2p^{n + 1)2 - 4/3(n + l)f'^ + t'* - 4:(3^]e 
''8(3^{n + 1)3 + 8(3^{n + 1) = 0. 
For e << t we can then obtain 
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where 



r 3 /32 «2 -, -1/2 

J^(/3,n)=[l-4^(n + l)-4^ + 12^(n + l)2 

(33) 
For small fields (i.e. (3 « t'2), F{(3,n) « 1. Therefore, 
in the limit of low energies and small fields, the Landau 
levels should approximately depend on the magnetic field 
as B^^^, in agreement with the results of Ref. [9]. 



IV. NUMERICAL RESULTS 

Let us first consider the ABA case. Figure 2 shows the 
field dependence of the low- lying Landau levels calculated 
from Eq. (18), with Ui ^ U2 = U3 = 0, for n = 
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FIG. 2: (Color online). The lowest Landau levels as function 
of magnetic field for ABA-stacked graphene trilayers calcu- 
lated from Eq. (18), with Ui = U2 = U3 ^ 0, for n = 
(black solid lines), n — 1 (red dotted lines), n = 2 (green 
dashed lines), n = 3 (blue dot-dashed lines) and 4 (yellow 
dot-dot-dashed lines). 
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FIG. 3: (Color online). Energy spectrum as function of mag- 
netic field for ABA-stacked graphene trilayers, for Ui — 100 
meV, U2 = 50 meV and C/3 = 25 meV, n = (black dots), 
n — 1 (red squares), n — 2 (blue lozenges), n — 3 (green 
triangles) and 4 (yellow triangles). 



(black solid lines) , n — 1 (red dotted lines) , n = 2 (green 
dashed lines), n = 3 (blue dot-dashed lines) and 4 (yellow 
dot-dot-dashed lines). We find that for every value of 
n there are two different types of low-energy branches: 
1) those that depend linearly on the magnetic field (i.e. 
"bilayer-like" behavior), and 2) branches that display a 
B^/^ dependence ( "monolayer-like" branches). A third 
set of bilayer-like branches are found around E = ±ti. 



FIG. 4: (Color online). Low-lying Landau levels as function 
of the potential in the inner layer for ABA-stacked graphene 
trilayers, for n = 1 (red squares), 2 (green lozenges) and 3 
(blue squares) for B = 3 T, Ui = U3 ^ 50 meV. 
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FIG. 5: (Color online). Low-lying Landau levels as func- 
tion of the potential in the uppermost layer for ABA-stacked 
graphene trilayers, for n = 1 (red squares), 2 (green lozenges) 
and 3 (blue squares) for B = 3 T, U2 = 50 meV, U3 = 25 
meV. 



not shown in the figure. 

The effect of a potential difference between the layers 
on the energy spectrum as function of nragnetic field is 
shown in Fig. 3. As in the previous case, the figure shows 
branches corresponding to n = — 4. The potentials in 
the different graphene layers are C/i = 100 meV, U2 = 50 
meV and C/3 — 25 meV. As seen from Eq. (18), for S — > 
we have solutions corresponding to E = Uj, j = 1,2,3 
and E = {Ui + U3)/2. Thus, we find that the "monolayer- 
like" branches are shifted creating a gap with magnitude 
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FIG. 6: (Color online). Landau level spectrum for the trilayer 
graphene for the ABC stacking, as function of magnetic field, 
with Ui = U2 = U3 = (a), and Ui = U3 = 50 meV, U2 = 100 
meV (b) for n = (black solid lines), n = 1 (red dotted lines), 
n = 2 (green dashed lines), n = 3 (blue dot-dashed lines) and 
4 (yellow dot-dot-dashed lines). 
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FIG. 7: (Color online). Landau level spectrum for the trilayer 
graphene for the ABC stacking, as function of magnetic field, 
with Ui = 100 meV, U2 = 50 meV, and U3 = 25 meV, for 
n — (black solid lines), n — 1 (red dotted lines), n — 2 
(green dashed lines), n = 3 (blue dot-dashed lines) and 4 
(yellow dot-dot-dashed lines). 



U1—U3, whereas for the "bilayer-like" states a smaller gap 
opens with magnitude {Ui + U3)/2—U2- One consequence 
of this difference is the appearance of level crossings as 
the magnetic field is increased. 

Figure 4 shows the low-energy Landau levels as func- 
tion of the potential in the inner layer, for the ABA case. 
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FIG. 8: (Color online). Landau level spectrum for ABC- 
stacked trilayer graphene as function of Ui, for B = 3 T and 
U2 = U3 = 50 meV, with n — (black dots), n = 1 (red 
squares), n — 2 (blue lozenges), n = 3 (green triangles) and 4 
(yellow triangles). 
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FIG. 9: (Color online). Landau level spectrum for ABC- 
stacked trilayer graphene as function of U2, for B = 3 T and 
Ui = U3 — 50 meV, with n — (black dots), n = 1 (red 
squares), n — 2 (blue lozenges), n = 3 (green triangles) and 4 
(yellow triangles). 



for n — 1 (red squares), 2 (green lozenges) and 3 (blue 
squares) for B = 3 T, C/i = C/3 = 50 meV. Notice that: 1) 
the lowest energy levels depend linearly on U2 for small 
U2 (<< Ui = C/3) and for large C/2 (>> Ui = C/3) values, 
and exhibit an anticrossing behavior for U2 ~ Ui — U3] 
2) the higher energy states are very weakly affected by 
the bias. A different behaviour is observed when one 
varies the potential at the uppermost layer (C/i), as seen 



in Fig. 5. In contrast with the previous case, the bias 
is seen to cause a significant shift also on higher-energy 
Landau levels. 

The energy spectrum for the ABC case is shown in 
Fig. 6a, as function of magnetic field and with Ui — 
U2 = Us = 0, ioT n = — 4. For small B-values we 
see a doubly-degenerate branch with E = 0, and a B^l"^- 
behavior for the remaining branches, which turns into a 
linear behavior at large E. In comparison with the previ- 
ous case, the results in the ABC case show the presence 
of pairs of branches at low energies, whereas in the ABA 
case one finds two sets of energy levels for each Landau 
index. That is caused by the fact that, in the ABC case, 
the remaining four branches are found around E = ±t, 
with t ss 400 meV. 

Figure 6b shows results for an ABC TLG with Ui — 
Us — 50 meV, whereas U2 — 100 mcV. In this case, the 
main effect of the potential difference is the lifting of the 
degeneracy of the n = state and a shift of the whole 
spectrum to lower energy with increasing magnetic field. 

In contrast. Fig. 7 shows the LL spectrum for Ui = 
100, U2 = 50 and U3 — 25 meV. The inset shows an 
enlargement of the region around E — 100 meV. In this 
case, the bias creates an energy gap, which can be found 
by setting /3 — in Eq. (30), which leads to solutions 
with E = Ui and E — U2- Notice also the existence of 
level crossings, as well as the peculiar small magnetic field 
behavior where there is a reversal of the ordering of the 
Landau levels as compared to the regular high magnetic 
field behavior. 

Results for the dependence of the energy spectrum on 
Ui is shown in Fig. 8, for S = 3 T and U2 = U3 = 50 
meV. As seen, the degeneracy of the n = is lifted for 
Ui ^ U2,Us- Moreover, when the magnitude of the po- 
tential in the uppermost layer is increased, the Landau 
levels tend to become degenerate. A quite distinct pic- 
ture emerges if one varies instead the potential in the 
middle layer (C/2), as shown in Fig. 10, for i? = 3 T and 
Ui — U3 — 50 meV. In contrast with the previous results, 
the spectrum shows a linear dependence on the potential 
and there are no degeneracies for the different Landau 
indices. As in the previous figure, a single Landau level 
at ii^ = 50 meV is found to be unaffected by the bias. 



V. CONCLUSIONS 

In summary, we obtained exact analytical expres- 
sions for the Landau level spectra of trilayer graphene, 
within a model that took into account the layer asymme- 
try induced by different electrostatic potentials in each 
layer. The expressions were obtained for both the Bernal 
(ABA) and rhombohedral (ABC) stackings, which were 
found to display quite distinct behaviors. As shown in 
previous work, the Landau level spectrum for the ABA 
case in the absence of electrostatic bias between the layers 
shows both a monolayer-like as well as bilayer-like charac- 
ter, indicated by the different magnetic field dependence 
of the spectrum. The addition of a potential difference 
between the layers shifts the spectrum and creates a tun- 
able gap between the electron and hole states, the size of 
this energy gap being different for the "monolayer" and 
the "bilayer" energy levels. Level crossings between the 
"monolayer" and "bilayer" Landau levels are found for 
certain values of the magnetic field. 

For the ABC case, the Landau levels have a magnetic 
field dependence which, in the absence of bias, has a 
B"^' ^-dependence for low energies. For stronger magnetic 
fields the Landau levels exhibit a linear i?-dependence. 
The introduction of electrostatic bias in the system lifts 
the degeneracy of the n = levels and creates a tun- 
able gap. The results show also the existence of level 
crossings at small magnetic fields. This model can be 
refined by taking into account second-nearest neighbor 
terms, as well as remote coupling between the lowest and 
uppermost layers. However, these additional terms are 
expected not to infiucncc the qualitative agreement be- 
havior of the present results. 
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